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We evaluate the Casimir force between two parallel plates in Randall Sundrum (RS) scenarios 
extended by q compact dimensions. After giving exact expressions for one extra compact dimension 
(6D RS model), we generalize to an arbitrary number of compact dimensions. We present the 
complete calculation for both two brane scenario (RSI model) and one brane scenario (RSII model) 
using the method of summing over the modes. We investigate the effects of extra dimensions on 
the magnitude and sign of the force, and comment on limits for the size and number of the extra 
- - - dimensions. 
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^—i ■ I. INTRODUCTION 

(n: 

' Theories with extra space dimensions have a long history in physics, at the beginning being motivated by the 

r^ , attempt to unify the fundamental interactions, as in the Kaluza-Klein theories [l|. Recently, they have enjoyed a 

fS i' resurgence with the development of models with warped extra dimensions [3] , which provide an explanation of the 

JL I large gap between the Planck and the electroweak scale (the hierarchy problem) and the cosmological constant. The 

Q\ • possibility of the existence of extra dimensions, as well as their geometry and size influence the structure of the 

r^ ' vacuum, in particular the evaluation of the vacuum zero-point energy, known as the Casimir effect [3J . The research 

'— ' . in this area is motivated in two directions. First, developments in the fundamental area of the structure of the vacuum 

^^ ■ quantum field theories have been extensively explored, with a view to understand the implications of extra dimensions 

►^ . [j|. Second, several measurements of the attractive force between parallel plates (and other geometries) have firmly 

(^ ' established the existence of quantum fluctuations. The level of precision reached by these experiments may be sufficient 

■^ , to test models with different geometries p . Casimir forces in space-time with D > 4 Euclidean dimensions were first 

■TJ" ' calculated systematically by Ambj0rn and Wolfram [0]. Previous studies have analyzed cosmological aspects of the 

^^ [ vacuum, such as the cosmological constant as a manifestation of the Casimir energy during the primordial cosmic 

f~^ ■ infiation [7|. The Casimir energy has been investigated in the context of string theories [8|, and even in Randall 

f^ ' Sundrum models, as means of stabilizing the radion 0- The dynamical Casimir effect has also been discussed in 

OO r warped braneworlds [Ifll • The Casimir force for parallel plate geometry has been calculated in UED [III, [Ij, [3 ^^^ 

>— ' ' in various other frameworks [l4| . 

The presence of the extra dimensions modifies the dispersion relation — r- ~ k\p, 4- Afc^p, where k^o is the usual 

. , 4 dimensional wave-vector and Afc^j-, would depend on specifics of the model used. By applying strict experimental 
jrt ' constraints to the resulting energy density or force, one can restrict the size or curvature of the new space dimensions. 
In a previous work [1^ , two of us analyzed the effects of one extra dimension on the Casimir force and energy in 5D 
dimensional warped space time of models RSI and RSII. In this work we extend our 5D calculation to Randall-Sundrum 
models extended by one, and then an arbitrary number of extra compact dimensions^. Such higher dimensional 
scenarios [16[ are interesting because they allow for localization of gauge and matter fields. For the extra spatial 
dimensions a common choice is to work in polar coordinates {Rg,6), where < i?6» < cxd and < 9 < 2tt. We expand 
and generalize the calculation of the Casimir force in the Randall-Sundrum model with two branes (RSI) and and 
one brane (RSII) between two parallel plates with arbitrary numbers of compact dimensions. We use as before the 
summation of the modes method. As far as we know, there are no alternative investigations of the qD RSI model in 
the literature. For the one-brane case (RSII) there exists an alternative evaluation using the Green function formalism 
[l7j . We compare our results with theirs and comment on any discrepancies. 
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Additional dimensions in 5D RS models have to be compact. Otherwise one could not recover Einstein's equations. 



Our paper is organized as follows. In the next section, we introduce some general features of the model in 6 
dimensions and derive the dispersion relation. In section III we present our analysis of the Casimir force in RSI 
models. We start with the 6 dimensional model (one extra compact dimension) and then generalize to an arbitrary 
number. We do the same for the RSII model in section IV. The numerical and graphical analysis of our results is 
included in section V, where we try to restrict the size and number of extra dimensions. We compare our results with 
others, summarize and conclude in section VI. 

II. RANDALL SUNDRUM MODEL WITH EXTRA COMPACT DIMENSIONS 

The metric of 5D space-time Randall- Sundrum (RS) model with 4D Poincarc invariancc is 

dsjo = e-2i^l^l77^,dx^dx'' - dy^ , (1) 

where k^^ is the radius of 5D Anti-de Sitter AdSs space, y is the spatial orbifolded extra dimension and yy^j^ is the 
usual 4D metric with the convention diag {+, — , — , — ). 

The generalization to more than one extra spatial dimension can also be done, and as mentioned earlier, these 
dimensions must be compact. Localization of graviton on the brane motivates us to consider additional dimensions 
both compact and warped [l^]- The metric in generahzed (3 + g + 2) dimensions (where 3-1-1 are the usual space-time 
dimensions, there is one extra y spatial dimension, and the extra q are compact and warped dimensions) is [16| 

dslo+n = e^'^^l^l I Vt^udx^dx" - ^ R^dO^ \ - dy\ (2) 

where y is non-compact and 0j 's are compact extra coordinates running from to 2t: Rj . This metric can be a solution 
to Einstein equations in (5 + q) dimensions with negative bulk cosmological constant [0, [3 • I^^ ^^^^ 6D case, for 
n = 1 becomes 

dsjjj = e-2^l^l [ri^^dx^'dx'' ~ R^dO^) ~ dy^ , (3) 

where we simply use R and 6 for the extra compact warped dimension. Here the 6D metric qmn is 

dslo = gMNdx'^'dx'' , (4) 

such that gMN = e^^*^'^! (im(7(r;^,y, — 1, — e^^l^l)^. The Greek indices ^,v run over to 3 while Roman indices M,N 
run over to 5. 

As in our 5D calculation ^^ , we use the scalar field analogy to represent the photon in 6D. The field equation for 
a scalar field (f> in 6D with mass nig can be obtained from the Lagrange density C{(f>, 9a/0), 

£(0, 9m0) = V^ (lg'"'dM<l>dNC^ - ^mecj^A , (5) 

where g is the determinant of the metric gMN which is y-dependent. (This has discussed previously in [l9[ and we 
include it here for completeness.) The field equation for (j> becomes 

(e^'^l^lTy^'-S^a, - i?-2e2i^lyl92 _ ^/^-ig,, {^^dy) + mf^ c^ ^ 0. (6) 

By using the separation of variables for (f>{x'^, 9, y) — X{x'^)Y{y)<d{9), one gets for j/ > 

d^Yiy) ^^ dY{y) 



dy dy 

-I 
1 



b-k^^^+m'^e^^yY{y)-mlY{y) = {), (7) 



p.^^_^ 2 _^(^)^0 (8) 

' dx^'dx'' ^ii "^ > ■• "^ > 



e„((?) = -==£-'', n- 0,1, 2... (9) 

V 27r_K 



Note that g^'^ = e^^^v^ diag{'qt^'' ,-l,-e~'^^\y\). 



where rn'^^ r = m^ + n^ / R^ is the effective mass in 4D. Since we use the scalar field-photon analogy, we consider the 
massless scalar field''. 

Equation (2.7) can be solved directly using the master formula for the equation of the form ,20j : 

y" + ay' + (6e^" + c)y = 0. (10) 

The solution for (2.7) with wlq ^ 0, with an appropriate choice of parameters, becomes: 

Y = eJ"'y''^[CiJ,{-me''y) + C2Y,{-me''y)] (11) 

k k 



where J^(z) and y^(z) are the Bessel functions of the first and second kind and v = ^ ■\/25k^ + 4m|. For the massless 
case V = 5/2. 

We apply the Neumann boundary conditions dyY{Q) = = dyY{TTRQ), to obtain the masses m as solutions to the 
transcendental equation 



J3 (^)r3 (^e-^«) - y3 (^)J3 (^e-^0) = (12) 



Expressing the Bessel functions J„_|_i (x) in terms of elementary functions we get, in the limit K,m > and 7rki?o 3> 1 

tan -^ '— = -^ ^ ' , 13) 

This is the exact solution. Following the same technique as in the 5D case, for 7rki?o ^ 1 we can neglect Y3/2 for 
large arguments and reduce equation (J12p to 

J^(!^e"'^«'') = 0. (14) 



The asymptotic form of J^(x) = \/ -^^ cos(x ) is valid for any v as long as x ^ li^^ — 1/41 is satisfied. Then, 

V na;^ 2 4 

the basic idea is to set the argument of cosine function (2A^ — l)7r/2, A^ = 1, 2, 3..., so that cosine vanishes. Thus, for 
arbitrary v we have 

x-'-^-^-^i2N-lf- (15) 



and get x as 



x=[N+^^)n. (16) 



This reduces to x = {N + l/4)7r for z^ = 1 in 5D case. In 6D (or any higher dimensions), we obtain f = 3/2 which 
instead gives x — {N + l/2)7r for the zeros of (2.14). (See the Numerical Analysis section for further details.) Putting 
the argument of J3/2 from above and using index N for m, which is now discrete, we get 

ruN = {N + ^)TT]!ie-''^"'\ iV= 1,2,3... (17) 

This is the form of m which should enter the energy density equation. Note that ttin ~ for A^ = , which will be 
important later. Also, for the 5D case m^v ^ {N + -)7rke~'"'^« , A^ == 1, 2, 3... [ll|. 

III. THE CASIMIR FORCE IN THE RSI MODEL 

In RSI, there are two branes, localized at y = and y = L, with Z2 symmetry y <-^ ~y, L + y <-^ L — y. The 3-branes 
have equal opposite tensions. The positive tension brane has a fundamental scale Mrs and is hidden; SM fields are 



There are some issues with the massive case in 6D for RSII since there are no zero modes. See Il9l for details. 



located on the negative tension brane, which is visible. The exponential warping factor gives rise to an effective scale 
on the visible brane located at y = L 

Ml, = Ms [1 _ e-2kL] (^8) 

A mechanism is needed to recover 4-dimensional General Relativity at low-energy, and this corresponds to introducing 
an extra degree of freedom known as the radion. A suitable choice of L (often taken to be related to the compactified 
radius, L — ttRq) and k allows the Kaluza-Klein (KK) spectrum to be discrete, and the lowest masses to be of ©(TeV), 
which predicts different collider signatures at low energies from those of large extra dimensions. As consistency with 
the low energy theory requires Mp\ ~ Mrs ^ k, there are no additional hierarchies in this model. The scale of 
the 4-dimensional physical phenomena space transverse to the S*'' dimension is then specified by the warp factor : 
A = AIpi e-^^'^o, and if we take A ~ 1 TeV, then we expect ki?o -- 12. 

A. The 6D RSI Case 

The effective dispersion relation in 4D can be obtained from (|5]) in fc-space. It becomes 

((^fc^)(^fc^)+m^^^)A(fc) = 0, (19) 

where k^ — (w/c, k). Hence, the following dispersion relation is obtained 
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«^(fc)=cEV''' + 152+"^e//: (20) 



If we further assume two-plates configuration separated by a distance a and apply Dirichlet boundary conditions along 
z-axis, we obtain the Casimir energy density per unit plate 



Ko - V.r.- j ^ l^pcj:^ i^l + — + ir.+ -^ I ' (21) 

where Vort is the volume of the orbifold and p is the polarization factor in 6D. This is our starting point in calculating 
the Casimir force. 

After finding the discrete m values corresponding to the RSI model with the above approximation, we can write the 
total frequency as 



^=p E "'y^l + ^ + ^ + '^'N (22) 

/.ri,JV=0 



where prime indicates that / = n = A^ = is excluded and tojv is given as 

if A^ = 



"^^~ ^ K(A + i) UN ^1,2,... 

with K = TT'ke~^^^'" . Note that the radius R of the sixth dimension is in general different than Rq of the fifth dimension. 
Rewriting the triple sum and rearranging the single and double sums to start with 0, we obtain 



e.n,N=0 



TT^P 



n 



2 



WlnN = pC ^ y ^'i + ~:2~ + "52 + ^N 



-^ J I L TTT^ 

i?2 



V- /,. 7r2^2 „2 V- /,, 7r2^2 „2 



i?2 

;,n=0 i,n,N=0 



N 



E/,T 7r2^2 ^2 ^2 

y'^^ + ^ + W^+T- (23) 



a- i?2 

!,n=0 



The Casimir energy density per unit plate is given by 

h f d^k 



^6D - ^orbf 



/i^rloV"^ 



7r2p „2 



a2 i?2 



"TV 



(24) 



We obtain 






E 



2' ^e,n=0 
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£,n,Ar=0 



In order to evaluate this expression we introduce the Epstein Zeta Function Zjv(ai, . . . , ajv; s) [2l| : 

ZAr(ai .. .,aAr;s) = ^ (aiTli + . . . flArn^) ' 

rt-i ,....7lN — — QiO 

and the inhoniogeneous Epstein Hurwitz Function E^ [2l|: 



(25) 



i?^(s;ai, . . . ,aAr;ci,. . . ,cjv) = ^ 



ni ,...,n w— 



TV 



^aj(nj+Cj)' 



We can write the energy density between the plates in a compact form as 



E, 



'6D ~^ ^ orb 



Vo, 



87rr(-i) 



ff E 



U,n=0 



3 7r2 1 



3 7r2 1 



a2 i^j +^^(-i^^'i^'-'^0'0'2)-^^^(-i^^'i^^«'«) 



pc;ir(-f)/i 1 TT 



1 ,7r 



1 . 1 



3 7r2 J_ 



^-''^^^^TX^(i^2(;^,-;-3) + -Zi(-;-3) + -Zi(-;-3) + ii;3(--;^^ 

i^2-(-|S,^;0,o) 



(26) 



(27) 



where the last form is obtained computing the explicit analytic expansion for the sums. Using the reflection formula 
for the Epstein Zeta function Q 



, n^ r(H^) ^ ^ 1 



1 



Zp{ai,...,ap;s) = «! ^ ■ ■ ■ S ^r^TTI wm ' ^p(r"' ■ ■ ■ ' ~;P ~ ^)' 



7rP/2 r(f) ''^ai 



we may write the first double sum on the right side 



^2(-^,-;-3) + Zi(^;-3) + Zi(-;-3) 
Ha ri a 



2r(-|)7r7/2 



3Ra / a \ I , , n'^ , ' 
^Z.(i?,-;5)+^C(4) + ^C(4) 



where in the last term we introduced the Riemann Zeta Function [2l|: C(s) = X]^i "■ 
The triple sum E3{~^; ^ , ^ , k'^ ; 0, 0, ^) is evaluated using [21| 



i^3(s; 01,02, as; Cl,C2,C3) = Jy^ E 



(-ir 



m 



-aTCni-'^rn, Ci)T{s + m)E2{s + m; 02, 03; C2, C3) 



m=0 

I A^r(s-i) 



1 



o\/ T./ x-^^2(s- ■^;a2,a3;c2,C3) 

2 y ai 1 (s) 2 

00 00 



/ TT C0S(27rCl) v-^ V^ / , « 3 r TT nf v-^ / n9i 

' ^ ^ ni— In2,n3— ^ j— 2 



(28) 



to obtain 



3 7r2 J_ 



,„,. . „ 1, ,,,31,1, la r(-2) , x n 



2'i?2' ' '2' 2V^r(-|) 



S2(-2;4'«';0,i) 



i?2' 



The double sum E^. 



V4 



'2 » 2'a2'i?2^ 
Zeta function and its expansion |2l| 



a0Fr(-|) 
3 7r2 1 



E E ^-^ (^ + '^^(^ + ^)^) ^2 (2a^y5 + .^(A^ + \y) . 



; 0, can be regarded as a special case of the inhomogencous Epstein-Hurwitz 



T{s) 



(-l)"T(s + m) fa 

.02 



„^ , s a, ■^-^ I — i r"i I S -h 777.) / ai \ , ^„rlnn, 

E^{s;ai,a2;ci,C2) = ffr E ^^TT^ ^ ( 7^ ) Cff (-2m, ci)£;i/"^(s + m; 1; ca) 



1/2-s ^ p/ _ j_\ 

-^1 (s-l/2;l;c2) 



ai T{s) 



2n' 



oo oo 



cos{2'KCi)ai 



-s/2-1/4 



^ ^nr'/'NK + C2)2 + c]-/2+V4 



m— 1 712—0 



X K, 



f^TTTll I -3 

.1/2 --^\^a2{n2 + C2)2 + c 

( J J \ LXJ -j 

where K^(x) is the modified Bessel function K^(x) = ^ —, with /^(x) = 7^ —r- — 

finally obtained as 

V4 



^^^«!(. + .)!V2 



2s+i/ 



(29) 



The sum is 



-b; 



i J,_L;0,0) . ,„(„.0)i5*(-|i„0) + l-^It|E*(-2;i,;0) 



00 oc 



■EE^- 
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Combining the three terms given above, we obtain, for the Casimir energy between the plates: 



(30) 



Jn-crj — Vorb 



1 



87rr(-i) 1 2r(-§)^5 



- — Z2 i?, -; 5 + -jC 4 + — C 4 



3 1 

2';r2' 



1 la r(-2) 



Cff(0,0)ii;2(-:^;^,K2;0,-) + -^^;A_^i?2(-2;^,«:2;0,^) 



1 

i?2^ 



00 00 



n,A''=0 



1, 



;;^pjy E E /-ni^ + «'(" + i)=j ^M ^"'V i^ + «'(" + 2'' 



1, 



C.(0,0)i.r^(-|i„0)4-^^i.r^(-2;4;0) 



i?2^ 



,^||-(S4)-(-;M)}- 



(31) 



Before proceeding further the energy without plates must be subtracted. Let's call the energy density without the 

np ■ 



plates e^p . The renormalized total energy shift can be given as 



^I (rcn) 
^6D 



E:,^" = E^nA - e^^A a. 



-6D^ 



'-np - 



(32) 



where A is the area of one plate. The energy density without plates ef^^ is given by 



,6D _ T^ ' 
triD ~ Vorb' 



/l^r\l 



iV = 



k +;^ + "^' 









V)i=l ri=OJV=0 

This yields, in terms of the previously defined functions 



E i/x' + ^ + E E v/"' + li + -5, - S 1^"^ + ^ + T 



n=0 



^2_ 

i?2 



(33) 



,6D 



16~ 



3 1 



C(5) + J:^i^.(-2;i„.-0,i) 



r(-2: 

7rir(- 



TT^i' (- 



2;^,0) 



Then the total energy in the 6D RSI is 



E. 



I (re 
6D 



= ^V, 



or6' 



87rr(- 



2r(-|)7r^ 



3 Ra ^ ,^ a ^. 1 ,,,, TT'',,,, 
- — Z2(i?, -; 5) + -^C(4) + ^C(4) 



aV^r(-l) 

2 



n 

i?2 



AaVorb fipc 3 1 



oo oo 



EE^-Mi. + T 




2 1 ' 



C(5) + (terms independent of a). 



(34) 
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We note that this expression contains no divergent terms and is finite, as expected, because the divergent terms in 
e^^ cancel exactly the ones in Eg^^ . The Casimir force is given by the derivative of the Casimir energy -Bg)^°" with 
respect to the plate separation a: 



F&D 



da 



(35) 



We use the the definition of Z2 to calculate derivatives. The force becomes, using ([M|) . 



A 


- "-^1 


1 
27r3 


- — Z2{R,-;b)~ 

a TT TT 




- l>\ 


^ 00 00 / ^ 




. 00 OC 

- -y y 

a ^—' ^—' 

£.Ar=in=a 


£-1 


($.«"(«. i)' 



15 -Ra^ 



E 

n,A''— — 00 
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a^ 



7l2 1 



K,{2ai\l^ + K^iN+lr 



(36) 



where we use the identity dzK^{z) ~ ~K^^i{z) K^{z). This expression is our final formula for the Casimir force 

density (force per unit area) for the 6 dimensional RSI model. Comparison with the expression for the force in 5D 
models [l5| shows that the terms in the square bracket arise as contributions from the extra compact dimension and 
are called UED-like terms, as in [ll|; while the last two terms given as Bessel functions are the contribution of the 
5D warped coordinates and will be called Bessel terms. In the limit i? — > 00: 
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'16^1 2^ 
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00 00 . 1 \ 1 1 
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(37) 



The first term is the Casimir force without the extra dimensions, and the two Bessel terms coincide with those in the 
5D consideration [l5|, except for the argument of the Bessel function. The difference in arguments arises from the 
location of the zeros of the Bessel functions in 5 and 6 dimensions, as shown in the previous section, and N + 1/2 
has to be replaced with iV + 1/4 as well, since the solution for tojv in 5D is different from the one in 6D. Note that 
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Thus the expressions (2.18) from [T^ and ([57)) are consistent with each 



B. The gD RSI case 

We continue by generalizing this expression to q extra compact dimensions. We proceed as before. We calculate 
first the energy density between the plates, then the energy density without the plates. We use this to evaluate the 
energy shift and then compute the force. 

The energy density between the plates for q compact dimensions (gD) is 
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(38) 



which give the energy density as 
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We obtain, after expressing everything in terms of Epstein Zcta and Epstein Hurwitz functions: 
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(40) 



where ^ indicates the sum over a permutation of the set g-terms of j-^, -^, . . . , -^j, X)o-i indicates the sum over 



a permutation of {q — l)-tcrms of{-^,-j|-,...,-^}. and so on. Note that X^i ^i i^ ^^'^ s^'^ °'^'^'" ^ permutation of one 
term of {f,-^,...,;^}. 

The energy density without the plate ej^p is given by 
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which gives: 
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TT^^" 
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The renormahzed energy shift is then given by 



E. 



I (ron) 
qD 



,A - e^„^Aa 



(42) 



(43) 



As before, the divergent terms in the vacuum energy without plates exactly cancel the divergent terms in the energy 
density between the plates. This yields, for the energy density, in a compact form 
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I (rcn) 
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,pchi 1 ^r(^) 
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-. oo 

^E^ 



£=1 Ai,n2,...,n<,=0 7V=l \i=l 

(terms independent of a) 



1 2 



E E(Ei + -^(^+|)^)^M2«^^Ei^ + -^(^+^)^ 



(44) 



where {^,Ri, ..., i?<j}(fe) represents all possible sets with the k element chosen in each case. So, {^, i?i, ..., i?q}(o) is an 
empty set. In the formula above we highlighted the finiteness of the energy density by using the reflection formula, 
see and ([28]) . allowing us to write: 



1 
'RI 



Zq^-,...,^^4) 



^ 2 l[|i?, Z,(i?i,...,i?,;g + 4) 
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^Zg_i({^^,...,^^}; 4) _ ^^+4p(_2)^ 



^ ^^''^Y.^Rl,---,Rq}Zq^l{{Rl,....Rqy,q + i). 



9-1 



For calculating the force, —i— 



^I (rc„) 
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E 



£,ni,n2,...,ng — — oo 



2*2 
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2„,2 



Rfni 



R^q 
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we obtain 



— Zg+i{-,Ki,...,Kg,q + A) _ -2 — 



E 



^^''"'-f + i?x 



.,ni ,712 ,...,ng — — CX3 



^g"g 



Further, for — — } —^ ' ' Z„({ — , i?i, . . . , R„}; q + 3), we note that only one term that will be independent 

of a. Thus, we may write it as 

TT- 2_^ —^ '—q — '■ Zq{{-,Ri, . . . ,Rq};q + 'i) = J> ^-n Zq{{-,Ri, . . . ,Rq};q + 3) 



9-1 



TT2 



9-1 ;,{«!,. ..,n<j} = -oo 

Then we obtain the general expression for the force in qD RSI, in compact form 
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^^^-^i^-b' 
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^ ^2 



.E|-^^(A'+^)^ 



(45) 



In differentiation we always assume that, other than - term, there is always at least one term chosen from the rest, 
{Ri, R2, ■■■, Rq}. We include the term with fc = 1 in the sum with the convention that {Ri, R2, ■■■, Rq}/o) ~ 1 ^nd 

As a check on our result, we verify for the case g = 1 it reproduces the 6D force found in the previous section. We 
have, for q = 1 
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(46) 



67r' 



Observing that Zi{—\A) = rC(4), the expression (I46p is identical to ([36]) . Thus the qD result coincides with the 

6D in the appropriate limit. 



IV. THE CASIMIR FORCE IN RSII MODELS 



In the Randall Sundrum model II (RSII), there is only one positive tension brane. It may be thought of as a limiting 
case of the previous model, where one brane is located at infinity (_Ro ~^ 00) • If there is a warped extra dimension. 
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large scales at the Planck branc are shifted at the TeV brane and the relationship between energy scales is given by 

M|i = ^^ . (47) 

Limits on Newton's law set lower bounds on the brane tension and the fundamental scale of RSII [22|. It has 
been shown also [23| that massless bulk scalars in the Randall Sundrum background (2.1) have similar properties 
as gravitons: there exists a localized zero mode and a Kaluza-Klein continuum of arbitrarily light states weakly 
interacting with matter residing on the brane. The spectrum of RSII is a continuous spectrum of to > KK modes, 
and there are no ©(TeV) signatures for this model at the colliders. The infinite dimension makes a finite contribution 
to the 5-dimensional volume because of the warp factor, and the effective size of the extra dimension probed is 1/k. 

A. The 6D RSII case 

We use the same procedure to evaluate the force in the RSII model as in the 5D case. The difference here is that 
the spectrum is continuous in to and consists of all tti > 0. The extra dimensional summation in N becomes an 
integral with measure dm/'k. We note that in RSII in addition to the continuous mode, there is the zero (discrete) 
mode which needs to be handled separately. The total energy density is the sum of zero and continuous modes: 



E^D,tot=E^D(o)+E^L- (48) 



where 



jii _ cp?i 

■"6£l(0) — ~^ 



,11 cpfi 




K'n - ^ I — I 7:^\ y.\IK + — + i^+m^\ (49) 

where the first expression is the discrete zero-mode and the second is the continuous to > mode. We deal with the 
zero-mode first. ^i^^/Q) is already calculated in RSI. It can be written as 



We re-write the zero mode energy density by using the reflection formula (j28[) to give 

We calculate the no-plate form of the zero-mode energy as well. The expression is the same as the zero-mass term in 
the no-plate energy density in the 6D RSI model: 



e 



6D(0) _ 



1 pch 3 2 



C(5). (52) 



"P 4 167r5 47r i?^ 

The energy shift for the zero mode becomes (including the area A and volume (aA) factors for the energy) 

pll (ren) q t \ 

^6D(o) pch (3Ra a c\ , ^ ._, , ^^ R^/-t.^ 3 2a \ 

^- = - 64^ U V^^(^' ^' ') + i?^^^ + ^^)^(') ~ 4^^ iF^^'V ^ ^ 

yielding the force for the zero mode component 

^ = 64^i47^^(^'^'^)-T— ^ nli^-^ + R^nlj - ^C(4) - ^^C(5) j. (54) 



ni,n2 — — 00 
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The continuos mode energy density (4.3) can further be written as 
cph f dm f cP]i.± 






{2ny 
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Ln=l 
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i?2 



n=l 



Ev/i-i 
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;=i 
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yielding, for the energy between the plates, 
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32 kV^r( 

The energy without the plate in RSII models is simply, for continuous modes 

cph f dm f d^h 
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i:\^^^+^ 



cph 



— 2C(6) 



i?2 

cph 



C(6). 



647ri?5ij(20F)7r¥"^'"' 327r7i?5k'^ 
We now compute the force in Casimir 6D including the zero mode contribution as 
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This is our final expression for the 6D force in RSII models, where 
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Comparing this expression with the corresponding one in 5D, we have, for i? — > oo 



F^^ 
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'647r5 






cpft 1 
64 kTT 



^C(5) 



^--lV2£^(^) 



(55) 



(56) 



(57) 



(58) 



where in the last expression we used FnoRS 



6ch 

'327r2a4 



C(4) 



n^ch 



240a 



J which is exactly what we obtained in (2.20) of 



B. The qD RSII case 



As we did in 6D, we quote the zero mode renormalized energy and the corresponding force from qD RSI results. 
We can get them from ((44]) and (|45)) . respectively. It is sufficient to switch off the term with Hurwitz Zeta functions. 
We have 
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II (ren] 
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pch 
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(terms independent of a). (59) 
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and 
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h 1 ''^^ p/^ fc+3 \ / 

= Oq+4 5 Yl 1~ Yl {Rl,---,Rq}{k-l) Zk{-,{Rl,...,Rq}{k-l);k + 3) 
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(60) 



Calculation of the energy and the force for the continuous mode gives 
cph f dm f (Pk. 



K'o-^-fl = 



pch 



For the energy without the plate in g'D for RSII models, we have 



oo q oc 
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(62) 



which yields 



.iD _ 
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Following the same method as used for the force in qD RSI model, we obtain the force 

A A da A da ' 

Then the general expression for the force in qD RSII models is, in compact form 
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(64) 



where k in the factor is the curvature of the space and the force corresponding to zero mode is (by reproducing the 
expression from (j60p ) 
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As a check on our result, we verify that in the case g = 1 it reproduces the RSII 6D force found in the previous 
section. We have, for g = 1 
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(66) 



Observing again that 2'i(-;4) = ^C(4), ^-^ 

TT a* 71-2 



c 2 



a2^2 



3 87r5 

4 a 



C(5)j El's well as 



RZi{R;6) — — ?C(6), the expression ([55)1 is identical to ([57]) . We have thus shown that the qD expression for the 
Casimir force reduces, with the appropriate substitutions, to the previously obtained expression in 6D. 



NUMERICAL ANALYSIS AND DISCUSSION 



In the previous two sections we presented complete compact expressions for the Casimir force between two parallel 
planes in Randall Sundrum Models extended by q compact dimensions. The models thus contain one non-compact 
extra dimension y, associated with a curvature parameter n which is determined by the 5D mass, as well as q compact 
ones. In the simplest case, the model will contain one compact (9) and one non-compact dimension (the 6D case). 
We analyze the 6D model first, then comment on an arbitrary number of compact dimensions by comparing 7D with 
6D. 

Additionally we distinguish between the two Randall Sundrum scenarios. In the (phenomenologically) more popular 
version of the model there are two branes, localized at y = and y — 7ri?o, with opposite tensions, called RSI. In the 
numerical analysis we concentrate on the RSI model, as RSII does not set any limits on the number or size of extra 
dimensions. The only approximation we have used in the derivation of the force is the form of rnjv. We obtained rriN 
as {N + 1/2)k in 6D by taking [N + 1/2)tt as the zeros of Bessel function J:^/^(xo). In 5D the relevant Bessel function 
was Ji{xo) and itln was approximated as {N + 1/4)k in our previous study [l5|. In Fig. [T] we show the difference 
between the exact zeros of Ji,3/2 and our approximation as a functions of N . As seen, the approximation works better 
as N increases, but is already good for small N values as well. 

The contribution of the compact dimensions (0i, 02, •■•) to the force resemble the one in universal extra dimensions 
(UED) , not surprisingly, as they enter into the dispersion relation in the same way. Whereas a full expression for the 
Casimir force for an arbitrary number of universal extra dimensions does not exist in the literature, our 6D result is 
the same as in [il| in the limit k ^ 0. We note that our force expression for qD naturally recovers the g-dimensional 
generalization of 5D UED as k —> 0. 

The difference between the Casimir force in RSI and RSII models rests in the appearance of the Bessel function 
terms, which are the contribution from the discrete KK spectrum, absent in RSII due to the lack of the second 
boundary condition. The Casimir force was calculated for RSII model in both 6D and qD case in [13| and in [rn |. 
respectively, using the Green's function method. We find it rather difficult to compare our results to theirs, as they 
do not present an expression for the force in terms of analytic functions. In the 6D analysis, they obtain the same 



P 
effective 4D Casimir force as we do (--FnoRs) 

both have F^^ 



under appropriate limits. Comparing their results in qD to ours, we 



FnoRs (1 + Aijs), but while in our case Ajis is a dimensionless term, as it should be, their term is 
of the form [(const)/Length], which is clearly problematic. 

We now proceed to analyze the implications of the expressions obtained on the size and number of the extra 
dimensions in RSI model. We first analyze the effect of one extra compact dimension in Fig. [2] as a function of plate 
separation a. We fix k = 10^^ GeV and the radius of 9 as Rg = 0.001 ^m in the left panel and Rg = 0.5 /xm in the 
right panel. In each case, there are curves for various kRy values. As we can see, all curves with ki?.y values less than 
21 coincide and are represented as one curve. In this case there is no contribution coming from RS (Bessel terms) 
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and our results reduce to those in to 5D UED. This feature can also be seen from our analytical expression given in 
([55)) . The RS effects are included only in the Bessel function terms and the rest is the same as for 5D UED for the 
6D case and for {q ~ 1)D for the qD RS case. We analyzed (p6|) more systematically using an analytical approach. 
This method is also helpful when estimating what one would expect numerically by adding more extra dimensions. 

We first tried to understand why effects appear from RS at around 'k.Ry = 21. We concentrate on the term with 
the modified Bessel function K2 (a similar analysis will hold for the Ki term as well): 



E^- 
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1, 
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K^{N 



1, 



(67) 



The dominant contributions come when the argument of K2 is small, thus the first few terms of series are dominant. 
In addition, as Re gets smaller and smaller, it increases the argument of K2 (for non-zero n of course,) and no matter 
how small k is, there is no contribution from terms with Bessel functions. Thus, the n = term in the series is 
dominant. For simplicity, take i = N = 1 and call the argument of K2 as Xn,e,N- Then, the dominant term is 
2^0,1,1 = 0,1^- The limiting form of Ka{x) is, for <C a; <C ya + 1 



K^{x) 



r(«) 



(68) 



The condition on xo.1,1 gives a critical value where the Bessel terms start being sizable. For example, for a — 0.5 /im 
and k = 10^^ GeV, the critical value of ki?j, becomes {'k.Ry)c = 21.0034. At any value greater than this critical value, 
K becomes exponentially smaller and the Bessel terms completely dominate contributions from the other terms. This 
is exactly the effect seen in Fig. [2] The same analysis for the Ki term gives 21.068 for the critical value of ki?j,. The 
limiting form of the Ki and K2 terms, including the factors in front (but not the sign) becomes [n/a)^ and 3/a^, 
respectively. 

We need to emphasize another feature. Like the Casimir force in 4D. the force in 6D vanishes as the plate separation 
increases. The C(5) term in ([55)1 is independent of a and blows up as Rg decreases. To see that the force in 6D behaves 
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FIG. 1; The approximations used for the zeros of the Bessel functions of first kind Ji{x) and J-i/2(x) for 5D and 6D, respectively. 
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FIG. 2: The Casimir force in QD as a function of plate separation a for various ki?y values with k — 10^^ GeV and Re = 0.001 urn 
in the left panel and Re — 0.5 /un in the right panel. The shaded region represents ±0.1 fim error for the measurement of plate 
separation a. 



like 1/a'* for very small Rg values, consider the Epstein Zeta function Z2{R,a/Tr;5) term in p6p . One can expand 
Z2{R, a/n; 5) as, defining /i = a/Rg, 



Z2{l,Vd;s) = 2as) + -^ 



20F r(^) 



d(s-i)/2 Y[s/2) 



C(.s-1) + Q(s,rf) 



where 



Q{s,d) 



, fjf , , , y n(-i)/2 y m- y cos f^^) A%-. (27rydn). 



In our case, d = (/^/tt)^ and s = 5. Thus, we have 
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Z2(l, m/tt; 5) = 2C(5) + ^C(4)— + Q(5, ^i'|^:') 



(69) 



(70) 



(71) 



If Rg is small for a fixed a, we have /i ^ 1 and ^'2(1, A'/"'; 5) goes to 2(^(5) which exactly cancels the other C(5) term 
in ([55)) . We arc left with the C(4) term which is very close to the 4D force term. Note that Q{s, d) also vanishes as /i 
increases. For example, the force in 6D becomes, as /i 3> 1 



Kd 



hcpn^ 1 3hcp 1 



{l + 2^^'^)Q{5,^lyn') 



(72) 



where the second term vanishes in this limit. This analysis can be generalized to the gD case and it can be showed 
that the same type of cancellations hold. Numerically, for small Rg values, we need the precise cancellation of very 
large numbers which require high numerical precision in the calculation. 

As can be seen from the above discussion, adding more compact dimensions would not affect much the kRy 
dependence of the force. Extra dimensions bring an extra suppression factor for the Bessel terms (when l/Rg^ and 
K are comparable) but these effects will be minor for most cases. The force will feel the effect of extra Rg's in qD 
coming from the terms other than the Bessel functions. As ki?j, is not affeted, in this limit the force is identical to 
that in UED in q dimensions. There exist some studies in the literature |ll| as well as some debate on whether the 
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FIG. 3: The Casimir force in 6D and 7D as a function of plate separation a for Rg — 0.001jj,m and 0.5 f^m with k = 10 GeV 
and k7?y = 15. Here Rg^ = Rg^ = Rg is assumed. The shaded region represents ±0.1 ^m error for the measurement of plate 
separation a. 

force is attractive or repulsive in extra universal dimensions [Ij] . As shown in previous sections, the 5D form of our 
force in K -^ limit is identical to the 5D result of [11| in UED and the force is always attractive. In Fig. [3] we 
display the force in both 6D and 7D as a function of a for two values of Rg, 0.001 fim and 0.5 fim. Here we assume 
Re^ = Re2 = Re- We chose ki?j, = 15, where there is practically no RS contribution from the Bessel terms so we can 
maximize the effects from extra compact dimensions. This can be considered like 6D UED in this limit. As in 5D 
UED, the force remains negative in 6D as well, thanks to the precise cancellation of a-independent 1/RJ terms (or Zi 
terms in general) with the ones from the Z2 term (similarly between Z3 and Z2 terms as well), as discussed above in 
6D RSI. In 5D UED, the bound from Casimir force on Rg is very weak (about 0.01/^m) [ll| but adding more compact 
dimensions, as seen from Fig. [31 will strengthen the bound. Conversely, if the limit on the size of extra dimensions is 
held constant, their number will be restricted. However, unless one introduces a large number of compact dimensions, 
it would be unrealistic from Casimir force considerations to get bounds on the size or number of extra dimensions as 
strong as the ones from b -^ 57, for example [2J|, which requires Rg to be less than 10^^^/im in 5D UED. 



VI. SUMMARY AND CONCLUSION 



We calculate, using dimensional regularization, the Casimir force between two parallel plates in Randall Sundrum 
models with extra compact dimensions. We include both the one brane scenario (referred to as RSII) and the two 
brane scenario (RSI). We present explicit analytic expressions for the energy between the plates, the renormalized 
energy, and finally for the Casimir force for the cases of one extra compact dimension (6D), and finally q extra compact 
dimensions (gD) in each model. 

We show consistency of the expressions both from qD to 6D and 5D models (RS models with no compact dimensions). 
We also show that in the limit in which the non-compact dimension of Randall Sundrum models becomes irrelevant. 
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our expressions reproduce exactly the known UED expression in one extra dimension. In that limit, our results present 
the generalization of the Casimir energy in q universal extra dimensions, and we settle the controversy regarding the 
sign of the force. We also compare our results with a calculation in RSII models in 6D and gD and briefly discuss the 
discrepancies. 

For the RSII model, comparison with the known Casimir force in 4D yields (as in 5D) negligible corrections, making 
this model not as interesting phenomenologically as RSI. In RSI, we perform an analytical and numerical analysis of 
our results and discuss the effects of extra compact dimensions on the size of ki?j,, the relevant parameter for RSI 
models, as well as the consequences of the RS models on the size of the compact (UED like) dimensions. Finally, 
we analyze the effects of having additional compact dimensions. Although the bounds on the size of the compact 
dimensions obtained are not as strong as those coming from other low energy experimental restrictions, our results 
include a comprehensive range of numerical and analytical explorations of the Casimir force in extra dimensions, and 
we are consistent with previous calculations in both RS and UED models in the relevant limits. 
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